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Abstract 

We present an approach for time-dependent quantum transport based on a self-consistent non- 
equilibrium Green function formalism. The technique is applied to a ballistic carbon nanotube 
transistor in the presence of a time harmonic signal at the gate. In the ON state the dynamic 
conductance exhibits plasmonic resonant peaks at terahertz frequencies. These vanish in the OFF 
state, and the dynamic conductance displays smooth oscillations, a signature of single particle 
quantum effects. We show that the nanotube kinetic inductance plays an essential role in the 
high-frequency behavior. 
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Nanoelectronic devices using nanotubes and nanowires as their active elements have been 
extensively studied for the, DC properties g However, their high-frequency AC character- 
istics have received little attention despite the obvious importance for many applications and 
the breadth of intriguing scientific questions. Experimental work in carbon nanotube (NT) 
and graphene field-effect transistors (FETs) have indicated little performance degradation 
up to GHz frequencies jsj; and recently, time-domain measurements in the terahertz (THz) 
regime have been used to distinguish between plasmon and single-particle excitations at low 
temperatures in NTs ^|. 

This progress in measuring the high-frequency properties of NTs poses challenging ques- 
tions for theory and modeling, in particular on how to describe the carrier quantum dynamics 
under non-equilibrium conditions in realistic device geometries, along with the self-consistent 
feedback between the time-dependent charge and potential, which is essential to capture the 
plasmonic excitations of the system J4]. In addition, one open question is whether single- 
particle and plasmon excitations can be distinguished in low- dimensional systems. 

In this paper, we address these questions by proposing a formalism for AC quantum 
transport making use of non-equilibrium Green functions (NEGF) [5], and apply it to de- 
termine the high-frequency properties of NTFETs. We find that the dynamic conductance 
exhibits both smooth oscillations and divergent features in the THz regime, indicating the 
coexistence of single-particle and collective excitations (plasmons). In addition, by calculat- 
ing the dynamic capacitance, we show that the nanotube kinetic inductance plays a central 
role in determining the high-frequency behavior, in contrast to conventional FETs. 

We begin by describing our theoretical approach using the example of a NTFET as 
illustrated in Fig. la. (NTs are ideal to study ballistic high-frequency transport due to 
the large electron mean-free path for acoustic phonon scattering, even at room temperature 
6[].) The system is divided into a device region which is connected to two semi- infinite leads 
consisting of NTs embedded in source and drain metals. The salient feature is the presence 
of a time-dependent potential at the gate terminal, V g + v g (t), and we are interested in 
the dynamic source-drain conductance g{uS). This is a much different problem than the 
one treated previously mostly within a non-self-consistent approach where either a time- 
dependent potential was app 



two-level systems were used 



ied to the source and drain electrodes [7| or simplified one- or 
8|- 



In the absence of v g (t), the system is described by the retarded/ advanced (7 = r/a) 
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FIG. 1: (color online) (a) Cross section of the cylindrical NTFET geometry, (b) DC transfer 
characteristics for a L = 30 nm long NTFET. (c) Band diagram in the ON and OFF states, as 
marked by arrows in panel (b). 

Green function 

Gl{E) = [(E ±irj)I-H - U sc - ^(E)]" 1 (1) 

where rj is a positive infinitesimal constant. H is the time- independent Hamiltonian of the 
device region while couples the device region and the semi-infinite leads; U sc is the time- 
independent, spatially-dependent, self-consistent DC electrostatic potential. This Green 
function forms the basis for calculations of the DC transport properties of nanosystems. 
In the presence of v g (t), a time- and spatially-dependent electrostatic potential U(r,t) is 
generated in the device region. U (r, t) must be determined based on the complex environ- 
ment of the nanosystem, and must satisfy self-consistency between charge and electrostatic 
potential, as will be discussed further below. Because of the presence of this time-dependent 
potential, the relevant Green function is the two-time function 

G 1 (t, t') = Gtf (i - + J dUGl (t - h) f/(t!)G 7 {t u t') . (2) 

We obtain the frequency-dependent particle current at terminal a from the time derivative 
of the number operator as 

JJ(M = \Tt I dE {G<(E + , E)K(E) - K(E + )G<(E + , E) + G r (E + , E)Z<(E) - E<(E + )G a (E + , E)} 

(3) 
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where E + = E + ftw. The AC Green functions G 7,< (£' + , E) appearing in Eq.Q are obtained 
from Dyson's equation G 1 = Gq + G^UG 1 , whereas the non-equilibrium particle density 
reads G< = (I + G r U)G£(I + UG a ) following the standard procedure fl. G< = G r ^<G a Q 
refers to the particle distribution of the unperturbed system with = — J2 a faO^a ~ ^<*)> 
and f a the Fermi function of terminal a. 

The above set of equations provides an approach to calculate the frequency-dependent 
quantum transport in the NEGF technique. These equations need to be augmented to 
include the coupling of the quantum transport equations with the electrostatics, as embodied 
by Poisson's equation: 



(At the frequencies of interest here, the time dependence of the full Maxwell equations can be 
neglected.) This contribution is especially important to capture the complex environments 
of nanoelectronic devices, and plasmonic effects. Poisson's equation is complemented by 
appropriate boundary conditions at the terminals; in particular at the gate terminal, a time- 
dependent potential v g (t) is applied and serves as the external perturbation that generates 



A closed set of equations can be obtained by expressing the charge density from the Green 
function: 



Thus, the set of Eqs. dH) — (El) provides an approach to calculate self-consistently the 
frequency-dependent current in the presence of a frequency-dependent external gate poten- 
tial. We note that our approach does not rely on approximations to the wide- or narrow-band 
limits, but treats the spectral properties of the contacts explicitly. 

To proceed further, we consider U (r, uj) to be a small perturbation and expand the above 
equations to linear order in U. The Green function of the perturbed system can be obtained 



from Dyson's equation as G 1 = Gl+^G2' + UG J , while G< = G^ + \G^ & UG%+\G r ^UG% (a 



+ superscript indicates a function evaluated at E + ftw). In practice, one is often interested 
in the small-signal response of the two-terminal conductance g a p (uj) = dl a (uj) /dVp\y Q . To 
obtain an expression for g a p, we consider a time-harmonic gate potential v g (t) = vqcos (ut) 
which generates an electrostatic potential U(r, uj) cos {ut) on the NT. Expanding to lowest 
order in U and Vp we obtain for the particle conductance 



V- [e(r)W(r,u;)] 



(4) 



U(r,u). 




(5) 
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<3 M = f^r | d£[{<^ + f/(u,)<^£|-£^ 

+G£ + C7 ( W ) G^a - T£+G$U (u) G a - Z r > + G r p ( w ) G^] (6) 

where fyg = d/i/cH^I^ for a general function h. 

In general, the particle conductance does not obey fundamental sum rules for current 
conservation and gauge invariance, since the displacement current has been omitted. To 
include it, we adopt the scheme of Wang et.al. [91], where the final conductance is given by 
9ap = 9ap- (Z] 7 = s ,dfi , Q 7 / 52s=B,d9s) 9% and derive for the frequency- dependent displacement 
conductance 

+G r ' + U(u)G r t<G«] . (7) 

Note that g a p{oj) measures the change in the conductance relative to the value at the DC 
operation point. 

We now apply the AC theory outlined above to calculate the zero-bias AC conductance 
at T = 300 K for the NTFET of Fig. la. The first step is to obtain the DC properties of the 



NTFET; details of the numerical procedures are given in Ref. [lOj. The (17,0) NT of radius 
R = 0.66 nm is modeled using a p 2 -tight-binding model with an overlap energy of 70 = 2.5 
eV, giving a bandgap of 0.55 eV. The Fermi levels Ep of the source and drain metals are 
set 1 eV below the NT midgap before self-consistency, which gives a p-type ohmic contact 
after self-consistency. Fig. lb shows the DC transfer characteristics of the NTFET with a 
channel length L = 30 nm, which consist of an ON-state where the bands are essentially 
flat, and an OFF-state where a gate-controlled barrier blocks the hole current (Fig. lc). To 
explore the AC behavior of the NTFET, we first choose a DC operating point, either in the 
ON or OFF state as marked by arrows in Fig. [lb, and apply an AC gate signal perturbation 
of frequency u and magnitude Vq = 10 meV. Figure [2^ shows the real and imaginary parts 
of the AC conductance in the OFF state which displays smooth oscillations as a function 
of frequency. (g(u = 0) is negative because the AC perturbation v cos (ut) reduces to a 
positive DC voltage perturbation 5V g = vq. According to the DC transfer characteristics 
of Fig. [TJd an increase in DC gate bias leads to a reduced conductance.) Surprisingly, for 
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FIG. 2: (color online). AC Response in the OFF-state for a L = 30 nm NTFET: (a) Real/imaginary 
part of g(uj). (b) Color plot of the density of states along the channel, showing the resonant 
photoexcitation of carriers through energy and spatially oscillating quantum states. The valence 
band edge is marked with the solid black line, and Ep = — 1 eV. 

finite frequencies the AC conductance can take values larger than the DC conductance. The 
origin of this behavior and of the smooth oscillations can be understood from the spatially 
and energy dependent density of states (DOS) in the OFF state shown in Fig. [2)3. At a given 
position along the NT the DOS shows oscillations in energy with a characteristic frequency 
of about 25 THz. The maxima in g(co) arise from the photoexcitation of carriers between 
maxima of the DOS while the minima in g(u) arise from the excitation between maxima 
and minima of the DOS. Thus, in the OFF state, smooth conductance oscillations are a 
signature of the single-particle excitation spectrum. The oscillatory character of g(u) is 
preserved when the self-consistent feedback between charge and potential is disabled (not 
shown), confirming this nature of transport. 

Figure 3a displays the real and imaginary parts of the dynamic conductance g(uj) for 
the same NTFET in the ON-state, which is also slightly negative at uj = for the same 
reason as in the OFF-state. We observe that the response for frequencies less than about 4 
THz is constant, in agreement with the experiments in Ref. Qj. For larger frequencies we 
observe that the conductance exhibits a pronounced divergent response at uj/2ti « 24 THz 
with an underlying oscillatory behavior. Near the divergence, the self-consistent potential 
along the NT reveals large amplitude oscillations, which change phase as the divergence is 
crossed. These large amplitude oscillations are shown for a device with L = 90 nm channel 
length in Fig. 3b, including the higher-order modes. The presence of these oscillatory modes 
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FIG. 3: (color online) Panel (a) shows the frequency-dependent conductance for a NTFET with 
L = 30 nm. Top inset compares the self-consistent and non-self-consistent results. Bottom inset 
shows the dielectric function versus frequency. Panel (b) shows the periodic potential near the 
plasmon frequencies for L = 90 nm. Panel (c) shows the response for three different channel 
lengths. 

suggests that plasmons are responsible for the divergent behavior. In order to ascertain that 
the response calculated in the ON-state can be attributed to collective rather than single 
particle excitations, we compare in Fig. [3^ (top inset) the self-consistent (SC) response 
g(oo) with the non-self-consistent solution obtained from the first iteration. The divergence 
in the conductance disappears entirely while the smooth oscillations persist. Thus, the 
self-consistency between the charge and potential is essential to observe the divergence, a 
signature of a collective phenomenon. In the ON state we therefore have the coexistence of 
single-particle and plasmonic effects. 

Furthermore, from the amplitude of the potential oscillations we can obtain the dielectric 
function of the NTFET from e{uj) = U ex t (uj) /U (u). As shown in the bottom inset of Fig.[3K, 
e{u) shows a clear zero crossing at the frequency where g (lo) displays divergent behavior, 
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and is further evidence for the excitation of plasmons. In the simplest model, e (q, u) = 
1 — x (q, K (q) where x (Qi w ) is the response function and K (q) is the electrostatic Green 
function for Poisson's equation. In the ON state, the NT is effectively like a quasi-one- 
dimensional electron gas, where for small q, xiQ^) < 2q 2 kF/miiu! 2 , with hp the Fermi 



wavevector and m the effective mass 



111 ]. K(q) is obtained by using the Green function for 



a NT in free space surrounded by a dielectric e ox , K{q) = j—Io{qR)K (qR) . (This is a good 
approximation to the full K(q) when q ^> Rq 1 where R g = 4 nm is the gate radius.) Near 
the threshold voltage we obtain the plasmon frequencies for the NTFET by quantizing the 
plasmon excitations to the channel length, i.e. q n = (vr/L) (2n + 1), giving the frequencies 
u n = -jj (2n + 1) I (nnR/L) KqIq (rnrR/L) where v p is the plasmon group velocity. To test 
the scaling of the plasmon frequencies with channel length and mode number, we calculated 
the AC response for NTFETs up to channel lengths of 90 nm; inspection of Fig. [3b indicates 
that the plasmon frequencies scale linearly with L~ l and n, with a plasmon velocity of 
v p = 2.7 x 10 6 m/s, about four times larger than the Fermi velocity. 

Our modeling approach also allows the study of the fundamental processes that control 
the AC properties of nanoelectronic devices. As an example, we show in Fig. H] the real 
part of the total dynamic capacitance defined as C{uS) = Q(oj)/vq with Q the total charge 
on the NT. At small frequencies Re[C(u;)] is positive implying a capacitive-like behavior, 
but becomes inductive at ~ 8 THz as marked by the sign change. This characteristics is 
fundamentally different from that of traditional FETs which show only capacitive behavior. 
The origin of this unconventional behavior is the nanotube kinetic inductance. Indeed, we 
can model the NTFET as a classical RLC circuit as shown in Fig. H] (inset), for which the 
dynamic capacitance Re[C(u;)] oc 1 —u 2 LkCq where is the nanotube kinetic inductance 
and Cq is the zero-frequency capacitance. This gives a transition to inductive behavior at 



oo = 1 /^LxCq, and from our numerical data for C{uS) we extract = 0.2 nH for the 30 
nm device. This value is consistent with that expected from simple arguments 12]. 

In summary, we presented a new self-consistent approach for AC quantum transport 
and applied it to determine the high-frequency response of NTFETs. In the ON-state, the 
dynamic conductance shows divergent peaks, which are associated with the excitation of 
plasmons of the gated NTFET acting as resonant quantum cavity, whose mode spectrum 
can be tuned by varying the channel length. Our results suggest that low-dimensional 
systems with nanometer sized channels show potential for novel detectors and emitters of 
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FIG. 4: (color online). Real part of the dynamic capacitance C(uj) for L = 30 nm in the ON state. 
Inset: RLC-model for the dynamic capacitance C(u>). 

terahertz radiation. The approach can be applied to a broad range of nanoelectronic systems; 
it will be useful for the study of many time-dependent phenomena in low-dimensionality 
systems including phonon and defect scattering, ultrafast optical excitation, and entirely 
new operation modes of nanoelectronic devices. 

We are indebted to M. Vaidyanathan and H. Guo for fruitful discussions. Sandia is a 
multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Co., for the 
United States Department of Energy under Contract No. DEAC01-94-AL85000. 
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